Structure Functions and Pair Correlations of the Quark-Gluon 

Plasma 



Markus H. Thoma 

M ax- Planck- Institut fur extraterrestrische Physik, P.O. Box 1312, 85741 Garching, Germany 



I. INTRODUCTION 



in 
o 

o . 

<N . Abstract 

> 
O 

^ | \ Recent experiments at RHIC and theoretical considerations indicate that the 

\@ ' quark-gluon plasma, present in the fireball of relativistic heavy-ion collisions, 

might be in a liquid phase. The liquid state can be identified by charac- 
■ teristic correlation and structure functions. Here definitions of the structure 

functions and pair correlations of the quark-gluon plasma are presented as 
\ well as perturbative results. These definitions might be useful for verifying 

the quark-gluon-plasma liquid in QCD lattice calculations. 
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Relativistic heavy-ion collision experiments at RHIC have found evidences for a new state 
of matter, the quark-gluon plasma (QGP) [1,2]. It is expected to exist in the hot phase of 
the fireball in ultrarelativistic nucleus-nucleus collisions as well as in the early Universe 
for the first few microseconds after the Big Bang at temperatures above T c = 170 MeV 
[3]. In the limit of an infinitely high temperature the QGP should be an ideal gas because 
the effective temperature dependent coupling constant becomes small due to asymptotic 
freedom. Therefore in the (extremely) high temperature regime the interaction between 
quarks and gluons is weak and the QGP can be described by perturbative methods [4]. 
However, at temperatures which can be realized in accelerator experiments, i.e., at maximum 
a few times of the transition temperature, the QGP is a strongly interacting many-body 
system. 

In classical, non-relativistic electromagnetic plasmas one distinguishes between weakly 
coupled and strongly coupled plasmas. For this purpose the Coulomb coupling parameter, 
T = Q 2 /(dT), where Q is the charge of the plasma particles, d the inter-particle distance, 
and T the plasma temperature, is considered [5]. In the case of an one-component plasma 
with a pure Coulomb interaction the Coulomb coupling parameter corresponds to the ratio 
of the interaction energy to the thermal energy per particle. Most plasmas in nature and 
in the laboratory are weakly coupled, i.e., r <C 1. If T is not much smaller than one, the 
plasma is called strongly coupled. These plasmas can be in a gas, liquid, or even solid phase 
[5]. For T > 0(1) the plasma behaves like a liquid and for T > 172 the plasma particles 
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are predicted to arrange in ordered structures, the plasma crystal [5]. The latter state was 
discovered in so-called complex or dusty plasmas [6]. 

In real plasmas the Coulomb interaction is modified to a Yukawa type interaction due to 
screening. Then the Coulomb coupling parameter, defined above, is no longer the ratio of 
the interaction to thermal energy. An additional parameter, the ratio of the inter-particle 
distance to the Debye screening length, k = d/Xo, becomes important. If k ^> 1 the plasma 
is weakly coupled. Now the transitions to the liquid and solid phases are considered within a 
phase diagram in the T-zt-plane (see e.g. Ref. [7]). For larger values of k a higher value of V 
is required to achieve a phase transition to a more ordered phase (gas-liquid or liquid-solid). 
However, as long as k is of the order of one or smaller, the phase transitions are not much 
shifted to higher values of Y [7]. In addition, it should be noted that in strongly coupled 
plasmas the screening potential contains a long-range power law contribution [5]. 

In the case of a QGP the coupling parameter is estimated to be V = 2Cg 2 / (47rdT) = 1.5 
- 5 [8,9]. Here C is the Casimir invariant (C = 4/3 for quarks and C = 3 for gluons), 
d ~ 0.5 fm the inter-particle distance, and T the temperature assumed to be about 200 
Mev corresponding to a strong coupling constant g ~ 2. The factor 2 in the numerator 
comes from taking into account the magnetic interaction in addition to the static electric 
(Coulomb) interaction, which are of the same magnitude in ultrarelativistic plasmas. 

As a matter of fact, the true coupling parameter of the QGP might be even larger because 
the potential might differ from a simple Coulomb or Yukawa potential corresponding to a 
one-gluon exchange. Indeed the cross section enhancement by about a factor of 80 due 
to non-perturbative effects discussed in Ref. [10] indicates a stronger interaction potential. 
The fact that the cross section enhancement by a factor 2-9 due to modifications of the 
scattering for T =1.5 - 5 [9] does not explain the observed cross section enhancement suggests 
that a realistic coupling parameter could be up to an order of magnitude larger. 

Following the investigations on non-relativistic plasmas, we proposed that the QGP in 
relativistic heavy-ion collisions is a liquid rather than a gas [9] which is supported by exper- 
imental results [2] as well as other theoretical considerations [11-14]. Among those are an 
observed strong elliptic flow, a strong jet quenching and a fast thermalization corresponding 
to large parton cross sections, and a small viscosity indicating the QGP as an ideal fluid. 

The presence of a liquid QGP expands the phase diagram of strongly interacting matter 
by a new phase transition from a QGP liquid to a gas at high temperatures. Also a new 
critical point in the phase diagram should show up, above which only a supercritical fluid 
exists. Assuming that this phase transition occurs at V ~ 1 - the distance parameter in the 
QGP, k — 1 - 3, is sufficiently small - and that the temperature dependent coupling constant 
decreases logarithmically with the temperature, the transition temperature is estimated to 
be of the order of a few T c . Such a temperature could be reached at LHC, where the 
liquid-gas transition might occur during the expansion of the fireball [15]. 

The phase diagram of strongly interacting matter can also be studied within lattice 
QCD. However, the phase transition from the liquid to the gas phase at high temperatures 
is probably difficult to observe in lattice calculations [16]. However, correlation functions, 
which can be investigated by lattice simulations, can also provide information on the state 
of the QGP. 
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II. CORRELATION FUNCTIONS IN LIQUIDS 



Quantitative investigations of liquids are possible by considering correlation functions 
[17]. In particular, the pair correlation or radial distribution functions in coordinate space 
as well as the dynamic and static structure functions in momentum space show characteristic 
properties within the different phases. For example, in the case of a liquid the pair correlation 
and its Fourier transform, the static structure function, exhibit a pronounced peak and one 
or two small and broad additional peaks. The first peak corresponds in coordinate space 
to the inter-particle distance, which is fixed in an incompressible liquid corresponding to a 
short-range order. In the case of a solid crystalline phase, where also a long-range order 
exists, a number of sharp peaks are observed. In the gas phase, where no order is present, 
the pair correlation function shows no clear structures. 

Here we want to consider the static structure function, which shows a similar peak 
structure as the pair correlation function [17]. Furthermore we will also discuss the dynamic 
structure function giving additional information about the system. First we will define 
these structure functions. These definitions might be in particular useful for confirming 
and investigating the liquid phase of the QGP in QCD lattice simulations. To demonstrate 
the structure functions and for comparison with future lattice calculations we will calculate 
them perturbatively. Of course, in this case we do not expect any liquid behavior, e.g., a 
single clear peak, as the perturbative regime corresponds to the gas phase of the QGP with 

r< i. 

The static density-density autocorrelation function is defined as [17,18] 

G(r) = l Jd"r'(p(r + r',t)p(r',t)), (1) 
where N is the total particle number and 

P M)=f>-r,(t)] (2) 
i=i 

the local density of point particles at the positions ri at time t. 

The density-density autocorrelation function is related to the pair correlation or radial 
distribution function, 

1 N 

^ r ) = ^(E^ + r i- r j)) ) (3) 

by 

G(r)=g(r) + 5(r). (4) 
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FIG. 1. Sketch of the static structure functions vs. momentum in the gas and liquid phase in 
arbitrary units [5]. 

The static structure function, denned by 

S(p) = ^ (P(P)M-P)) (5) 
with the Fourier transformed particle density, 

p(p) = Jd 3 rp(r)e-^, (6) 
is the Fourier transform of the density-density autocorrelation function 

S(p) = J d 3 r G(t) e~ ip r . (7) 

For uncorrelated particles the structure function is constant for p ^ [18]. The typical 
behavior of the static structure function in an interacting gas and in a liquid is shown in 
Fig.l according to Ref. [5]. There a one-component Yukawa system was assumed. Note 
that in a one-component Yukawa system there is no gas-liquid phase transition but only a 
supercritical fluid exists. After all the structure function shows a gas like behavior for small 
coupling parameters, whereas a peak structure (fluid like behavior) develops only at rather 
high values of T [5] . Only in the presence of repulsive and attractive forces at the same time 
a gas-liquid transition shows up, ending at a critical point. In the QGP, where a gas-liquid 
phase transition could exist [19], the appearance of peaks in the static structure function 
could correspond to a real liquid phase. 

The time-dependent generalization of the density- density autocorrelation function, 

G(r,t) = l |dV(p(r + r',t)p(r',0)>, (8) 
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is called van Hove function and its Fourier transform, 

S(p, u) = JdtJ d 3 r G(r, t) e ^'~ p - r) , (9) 

dynamic structure function. It is convenient to decompose the van Hove function into a self 
and a distinct part 

G(r,t) = G a (r,t) + G d (r,t), (10) 

with 

G.(r,t) = ^(i;«5[r + r I (0)-r I (t)]>, 

i=i 

G d (r,t) = ^(f:5[r + r i (0)-r j (t)]). (11) 

The distinct part looks similar as the pair correlation function, i.e., showing a pronounced 
peak plus one or two broad side peaks in the liquid case, if t <C r where r is the characteristic 
relaxation time of the system. For large t, however, it becomes a smooth and flat function 
of r. 



III. STRUCTURE FUNCTIONS IN A QUARK-GLUON-PLASMA 

Now we want to define these correlation functions in the case of the QGP. For simplicity 
we consider here only the quark component. The local quark density is given by 

p(r,t) =q(r,t) lo q(r,t), (12) 

where q and q are the quark wave functions. In a homogeneous and isotropic plasma, in 
which (p(r)p(r')) depends only on r — r', the static quark density-density autocorrelation 
function is given by 

G(r) = i(p(r,t)p(0,t)>. (13) 

n 

Here n = N/V = (p(r)) is the average particle density in a homogeneous system. At finite 
temperature its Fourier transform x follows from [20] 

G(r) = - A lim T jgj ^ X K,p), (14) 

where cu n = 2nnT are the discrete Matsubara frequencies. 
The corresponding spectral function is given by [20] 

ff(w,p) = - Imx(o;,p), (15) 

71 

where w is a real number now. The static structure function follows from the spectral 
function according to [20] 
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1 f°° UJ 

S(p) = - dua(uj,p) coth— , (16) 
n Jo 21 

which is a consequence of the fluctuation-dissipation theorem [18]. 

The function x( w ? p) can be expressed by the one-loop polarization diagram containing 
dressed quark propagators [20] 

00 r d 3 k 

X (uj, p) = -2N C T ]T / Tr[ 7o S F (k , k) 7o S f (uj -k ,p- k)] (17) 

n=-oo J \ An ) 

with the number of colors iV c and ko = 2nmT . This expression is related to the longitudinal 
part of the QCD polarization tensor containing only the quark loop, ITl = noo> by 

AN 

x(w,p) = — ^n L (w,p), (is) 

where g is the QCD coupling constant. 

The longitudinal part of the polarization tensor also determines the quark number sus- 
ceptibility Xq(T), which is proportional to the p = limit of the static structure function 
(16), Xq{T) = (n/2T) 5(0) [21]. In the case of an isotropic and homogeneous plasma the 
polarization tensor depends only on uj and p — |p|. 

According to (9) and (14) for t ^ 0, the dynamic structure function is given by 

S(u,p) = - X (co,p). (19) 

71/ 

Using (14) to (19) the density-density autocorrelation and structure functions can be 
derived in the strongly coupled phase of the QGP by using lattice QCD and investigated 
for signatures of the liquid phase as discussed above. In lattice calculations the correlation 
functions can be derived directly whereas the structure functions and spectral densities 
require a maximum entropy analysis [22] . For comparison, we calculate in the following the 
static and dynamic structure functions within perturbation theory, i.e., in a weakly coupled 
QGP, where no pronounced structures in these functions are expected. 



IV. HARD THERMAL LOOP APPROACH TO THE STATIC STRUCTURE 

FUNCTION 



Let us first consider the high-temperature or hard thermal loop (HTL) limit [23-25] . In 
this case \ is given by the high-temperature limit of the 1-loop polarization tensor with bare 
quark propagators, 



2N r Nf 



oo , uj + p 
— In 

2p uj — p 



(20) 



where N c is the number of colors, i.e., N c = 3 and Nf the number of light flavors in the 
QGP. For the interesting regime of small times, t <^ t, corresponding to uj p in (20), x 
behaves as 
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HTT / 2N c N f 9 p 2 
X RTL (lu » p) = T ^ (21) 

and the dynamic structure function (19) is a monotonically decreasing function of uj. In the 
case of a liquid this function also decreases with u but shows a peak structure in addition 

[17]- 

The spectral density, following from the imaginary part of the polarization tensor ac- 
cording to (15), within the HTL limit is given by 

a(co,p) = MlT^6(p 2 -co 2 ). (22) 
V 

The imaginary part of (20) corresponds to Landau damping. The static structure function 
according to (16) reads 

2NtT 3 

S(p) = (23) 

Here we have used consistently the high-temperature approximation, coth(cj/2T) = 2T/u, 
and N c = 3. The high-temperature limit corresponds to the classical limit, in which the 
polarization tensor, related to the dielectric function (see below) , can also be derived by the 
Vlasov equation [4]. Using the HTL polarization function together with coth(a;/2T) in (16) 
leads to the unphysical result of a linearly increasing static structure function for large p. 
The constant structure function (23), on the other hand, indicates an uncorrelated QGP 
corresponding to an ideal gas - the corresponding pair correlation function (4) is equal to 
zero - in the limit of infinitely high temperature. This result also agrees with the free quark 
number susceptibility XgCO = NfT 2 [21]. Let us note here that using the 1-loop polarization 
tensor beyond the HTL limit [26-28] in (16) also leads to an unphysical result, namely a 
monotonically decreasing S(p). This shows once more that the naive use of perturbation 
theory for gauge theories at finite temperature is inconsistent [25]. 

To go beyond the uncorrelated QGP, we resum the HTL polarization tensor (18) within 
a Dyson-Schwinger equation, 

oo 

ni = nf TL £(/^n? TL r, (24) 

n=0 

where D\ = 1/p 2 is the bare longitudinal gluon propagator in Coulomb gauge [4]. Using the 
relation between the longitudinal polarization tensor and the longitudinal dielectric function 
[4], holding in particular for the HTL limit, 

£ HT V , p) = 1 -Hl^>, (25) 

the resummed polarization tensor is given by 

^•^Grbr 1 )- (26) 
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The static structure function, following from (15), (16), (18), and (26) together with the 
high-temperature approximation, coth(cj/2T) = 2T/u, using a Kramers-Kronig relation is 
given by [18] 



4N C T 



S{p) = — c — p 2 



1 - Re ' 



e* TI V = 0,p) 



2N,V p> (27) 



9 i 2 ° 

n p z + m D 



where m? D = — n^ TL (o; = 0,p) = Nfg 2 T 2 /Q is the quark contribution to the classical Debye 
screening mass. The static structure function (27) starts at zero for p = and saturates at 
the uncorrelated structure function (23) for large p. Such a structure function corresponds to 
a Yukawa system in the gas phase (see Fig.l). Indeed, the pair correlation function following 
from the Fourier transform of S(p) — 1, 



NrT 3 



rn 



2 



g(r) = ^ e~ m °\ (28) 

Aim r 

reproduces the Yukawa potential. The radial distribution function, giving the probability 
of finding other quarks, in Ref. [5] is given by 1 + g(r). Using (28) the radial distribution 
function becomes unphysical, i.e., negative, for small r indicating the failure of the Vlasov 
approach at small inter-particle distances [18]. 

In a next step one could use HTL-resummed quark propagators and quark-gluon vertices 
in the polarization tensor as done for the quark number susceptibility [21]. However, to 
look for a strongly-coupled liquid phase within the structure functions and pair correlation 
functions, one has to adopt non-perturbative methods. At the classical level molecular 
dynamics, generalized to the relativistic case, could be useful [5]. In general, of course, QCD 
lattice simulations would be the ultimate choice. 



V. SUMMARY 

In summary, we have presented definitions of the pair correlation function, the density- 
density autocorrelation function, its time-dependent generalization (van Hove function) and 
their Fourier transforms, the static and dynamic structure functions in the case of a QGP. 
The latter are closely related to the longitudinal part of the polarization tensor. In liquids 
these correlation functions and their Fourier transforms show a specific peak structure. 
Hence they could be used to investigate the existence of a liquid phase in a strongly coupled 
QGP, if calculated within non-perturbative methods, in particular lattice QCD. Adopting 
the perturbative high-temperature limit of the polarization tensor and its resummation 
corresponding to the weak coupling limit, we have shown, as expected, that there is no 
indication of a liquid phase in the dynamic as well as static structure functions. These 
results can be used as a reference for non-perturbative calculations. 

Acknowledgment: I would like to thank W. Cassing, A. Peshier, F. Karsch, and W. 
Zajc for helpful discussions. 
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